#ifndef __FUN_H__
#define __FUN_H__

#include "hscuda_rasterdata.h"

template<typename DataType>
void printmat(const hs::raster::RasterView<DataType>& mat)
{
    for (int y = 0; y < mat.sizeY(); ++y) {
        for (int x = 0; x < mat.sizeX(); ++x) {
            printf("%6.0f", mat(x, y));
        }
        std::cout << std::endl;
    }
    std::cout << std::endl;
}

float reduceSum(const hs::raster::RasterData<float> & src_datas);

void reduceMaxIndex(const hs::raster::RasterData<float> & src_datas, float* p_max_val, int* p_max_id);

#endif
